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We consider extensions of the twin-trap Bose-Einstein condensate system 

of Javaneinen and Yoo [Phys. Rev. Lett., 76, 161-164 (1996)] to include 
pumping and output couplers. Such a system permits a continual outflow of 
two beams of atoms with a relative phase coherence maintained by the detec- 
tion process. We study this system for two forms of thermal pumping, both 
with and without the influence of inter- atomic collisions. We also examine 
the effects of pumping on the phenomenon of collapses and revivals of the 
relative phase between the condensates. 
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I. INTRODUCTION 



Both before and since the recent demonstrations of Bose-Einstein condensation 
in dilute alkah gases, the concept of the phase of a Bose-Einstein condensate (BEC) has 
attracted a deal of theoretical study. Traditionally, the existence of a phase is taken as a 
signature of spontaneous symmetry breaking and strictly it is only the relative phase of two 
BEC's that can be assigned a definite value. Many people have discussed the difficulties 
associated with the fact that in many cases, we may know the number of atoms present 
fairly precisely. The condensate is thus in a state that cannot possess a well-defined phase. 
Recently, several authors have modelled systems containing two BEC's, demonstrating 
that a relative phase can arise naturally, even when both condensates are initially in number 
states. Typically, atoms are permitted to leak out from both traps and are detected by 
some apparatus below. As it is unknown from which trap any individual atom comes, the 
distribution of positions at which atoms are detected shows interference fringes. At the 
same time, the quantum state of the two traps evolves from a simple product of number 
states into an entangled state of varying number differences between the traps allowing a 
well-defined relative phase to appear between the two condensates. The value of the relative 
phase is randomly distributed from run to run so that for an ensemble of runs, the phase 
symmetry is restored. 

The infiuence of inter-atomic collisions on such a system has also been consid- 
ered [[7|,^, [lO| , |lT| . The notable results are a reduction in the visibility of the observed in- 
terference pattern, and in the conditional visibility of the entangled state, due to diffusion 
of the condensate phase. Moreover, stopping the detection process after a relative phase 
has been established leads to collapses and revivals in the conditional visibility [p|, p!0| , p!T| , as 
the collisions cause the phases of different components of the entangled state to rotate at 
different rates. As the total state is a discrete sum of number difference states, the phases 
realign periodically and the visibility is restored. 

An obvious consequence of detecting atoms is that the system can not attain a steady 
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state. In time, the trap occupancies fall so low that the entangled state is reduced in size until 
eventually the atoms in one or both traps are exhausted altogether. If one envisaged using 
the two-trap system as a real "device" rather than merely an analogy for a single condensate, 
this might be a serious problem. A natural use would be as a kind of two-beam "atom laser" 
in which additional atoms would be tapped off from each trap through an output coupler 
into separate beams (see Fig. jT]), as in the case of two remarkable recent experiments at 
MIT [|3,|1]. The relative ph ase between the EEC's in the traps built up by the measurement 



process would then be reflected in the same well-defined relative phase between the two 
output beams, that could be exploited for some other interference experiment. While such 
a scheme could operate in a pulsed fashion where the traps were repeatedly filled with 
atoms, measured until exhausted and then the process repeated (as was the case for the 
experiments just mentioned), it would also be useful to have a continuous output. This 
would clearly require a continuous pumping of new atoms into the traps. We note of course, 
that interference measurements on the output beams would themselves act to produce a 
relative phase between the traps. It may be however, that the desired rate of measurements 
on the output beams is too low to stabilize the phase over long periods, whereas the detection 
rate directly on the traps may be as large as necessary. Here we assume that the rate at 
which subsequent measurements are made on the output beams is so low that they have a 
negligibly small feedback on the entangled state in the traps. (Note that the output rate 
may be relatively high as long as most of the leaked atoms are not actually detected in an 
interference experiment.) 

In this paper, we explore the effects of pumping a two-trap BEC system. We investigate 
what kind of steady states may be reached when pumping is included, and study the compet- 
ing effects of the collisions and measurements in such a steady state. We also allow output 
couplers on each trap as discussed above. We explore two types of pumping from thermal 
atom baths coupled to each trap. In the first, we allow two-way pumping where atoms are 
exchanged with the baths in both directions. Such a scheme was considered for a single-trap 



atom laser based on evaporative cooling by Wiseman et al jlj]. In the second model, atoms 



may enter the traps from a thermal bath, but the reverse process is forbidden. This kind 
of irreversible pumping scheme has been considered in an atom laser model by Holland et 



al |T^. We can then consider our whole system as a means of transferring non-condensate 
atoms in a thermal bath, into two coherent beams with the coherence established by the 
detection process. In this sense, the system might be considered a primitive two-beam atom 
laser. We emphasize though, that in this paper we include no line-narrowing element, a 



central component of true optical lasers [[lq -[18[|. The linewidth of the output beams would 
be at least that of the output couplers. 

The paper is structured as follows. In section |I| we describe our model in detail while 



in section 111 we discuss the nature of the entangled state more fully. Using Monte-Carlo 



wave function simulations, we consider pumping from a thermal bath of atoms in section |V 
examining the visibility and other parameters, when the pumping of the trap is either 
two-way or inwards-only . Finally in section 0, we turn to the phenomenon of collapses 
and revivals of the condensate phase and present some interesting effects associated with 
pumping. 



II. ELEMENTS OF THE MODEL 



We now set out the model in full detail. In all the work below we assume a system 
of two traps containing condensates with occupation numbers ni and n2 (see Fig. |^). We 
occasionally write rii to indicate either trap. Atoms may leak from either trap and be 
detected below at a mean rate establishing a relative phase. (We assume that an atom 
in either trap has the same probability of detection). The traps are pumped from two 
thermal reservoirs containing A'^i and N2 non- condensate atoms with rate coefficients Xi and 
X2 respectively. At different points in the paper, we assume that atoms may move either in 
both directions between the traps and reservoirs, or only into the traps, so that we define 
separate rates x™, a-^id Xi"* and X2^^ ■ the simplest systems, we would expect xT = xT^-i 
but the inclusion of some irreversible pumping process can prevent atoms escaping from the 
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trap into the baths giving = 0. One example would be to couple the thermal bath to 
an excited trap level |e). An atom in |e) can decay to the BEG mode \g) by spontaneous 
emission. If the medium is optically thin, the emitted photon is lost and excitation out of 
the ground state is impossible. We see in section |V| that this one-way pumping leads to 
quite different behavior to the two-way pumping. The physical validity of the two types of 
thermal pumping has been discussed at some length in Ref. . 

We also allow a separate leak from each trap into an empty mode at rates vi and z/2 to 
act as "output" beams. While we here leave the mechanism unspecified, we note that several 



techniques for creating an output coupler have been demonstrated by the MIT group |jT2[ in 
which rf signals or a bias in the trapping field are used to couple a portion of the condensate 
to untrapped spin states [|19|. Finally, the atoms in each trap experience collisions amongst 
themselves at a rate k. 

The master equation for the system may thus be written 



2-K 



^='-[p,H]+^ (0)] d<Pp + xT (Ni + 1)D [ai] p + xTN,D 

+XT (^2 + 1)D [aa] p + X2N2D Ul p + ViD [ai] p + V2D [02] p. 



P 



(1) 



where the Hamiltonian describing collisions amongst the atoms is [O] 

H 



K 
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a\ai 



(2) 



and is the annihilation operator for an atom in trap i. Also, for an arbitrary operator c, 
the superoperator D [c] p is defined by 



D[c\p = cpc^ — - {c^cp + pc^c) 



(3) 



The field operator ip (0) = ai + 026"*'^, where (j) = 2ttx, describes the detection of an atom 
at position x. Most of our results are obtained from Monte-Carlo wave-function simulations 
of Eq. (|1|) in which all leaks and additions of atoms to the traps are represented as quan- 
tum jumps, and the non-unitary evolution of the wave function is given by the effective 
Hamiltonian 



7 (alai + ala2) + [xT i^i + 1) + X^iVi + u,] a\ar (4) 
+xrAri + (A^2 + 1) + X^"iV2 + z/2] 4«2 + X^'iVsj . (5) 



When the chosen jump is a detection, the phase of the next detection is chosen randomly 
according to the conditional probability distribution 

P{<P) = {t\^H<P)^P{cP)\t), (6) 

where \t) is the state of the system immediately preceding the detection. It has been shown 
elsewhere that this may always be written in the form 

P(0) ocl + /5cos(0 + ^), (7) 

where the conditional visibility (3 and conditional phase 6 are determined by the previous 
history of the system. 

In order that the average number of atoms in the traps become constant in time for 
the thermal pumping schemes, we require a relation between the various coefficients in the 
master equation (|1|). Assuming equal detection and pumping rates for traps 1 and 2, the 
pumping rates must satisfy 



for two-way pumping with x*^ = = x' 



out . 



^ _ 7 {ni + ^2) + CTi (rii) + 0-2 (^^2) 
^ ~ m +N2 -{ni + n2) ' ^ ^ 

• for one-way pumping with = X™; ^ind = : 

^ " iVi((ni) + l) + iV2((n2) + l) ■ ^ ^ 

Note that for Ni ^ N2 ^ (rii) , the two-way pumping rate x"^ is larger than the one-way 
rate x~^ by a factor of approximately {rii -|- ^2) /2. In the one-way case, on average one atom 
is added for each atom detected or lost to an output coupler, while in the two-way case, on 
average all the trapped atoms are exchanged with the reservoirs for every loss by detection 
or output coupling. This difference has important consequences below. 
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III. QUANTUM STATE OF THE TRAPS 



Our main interest in this paper is to find the equihbrium state of the model just described, 
under a variety of conditions and explore the different influences of detection, pumping 
and collisions. For example, we naturally expect increasing collisions to reduce the phase 
coherence and drive the state towards a narrower number distribution. In preparation, 
however, we should first discuss the nature of the entangled state and our methods for 
characterizing it more fully. 

In earlier studies that consider only detection of the atoms the initial state is 

normally taken as the product state | )q = \ni,n2) with ni atoms in trap 1 and n2 atoms 
in trap 2. The unnormalized state following a single detection with phase is found by 
applying the field operator ip 

I )^ oc (ai + a2e~"^^) 722) = y/n^ \ni - 1, 722) + e~"''^ y/n^\ni,n2 - 1) . (10) 

By extension, after m detections the state has the form 

m 

\)m = ^'^k{m)\ni-m + k,n2-k), (11) 

A;=0 

where the are functions of the phases of all the detected atoms {0i, . . . , 0^} • If collisions 
are included, the state experiences unitary evolution under the Hamiltonian (0) in between 
detections and the coefficients Ck {m) are also functions of time. Here, with the inclusion of 
pumping the situation is similar, but as one can continue to detect atoms indefinitely, the 
entangled state can become very large (note of course that the number of detections m can 
now arbitrarily exceed the initial occupancy of the traps). It becomes more natural to drop 
the dependence on m and write the state at time t as 

p 

\t) ^Y^Ck [t) Mt) - k, n2{t) + k) . (12) 

k=—p 

This is an approximate relation because we truncate the sum at some cut-off p. This is 
particularly important numerically as the exact state can become prohibitively large for 
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calculations. Such a truncation is possible because the probability of all detections occurring 
from a single trap is small (assuming the initial trap occupancies are not wildly different) 
and hence the coefficients at the extreme ends of the entangled state are negligible. In our 
simulations, we drop terms for which \ck\ < 10"^^. We characterise the state in terms of the 
mean number of atoms in each trap, and the width of the number difference distribution 
with the natural definition 

o-n = {{{ni -riif) - {rii -^2)^)^^^ (13) 

Frequently we also wish to describe the phase distribution for which we use the width 



(T0 = (1 - \ {exp {i4)))ff^^ 



1 




(14) 



For a minimum uncertainty state, we have cr„cT0 = 1. Aside from the evolution of the 
conditional visibility, our main interest below is in the behavior of these two measures of the 
state. 



IV. STEADY STATES 

We now turn to finding the steady states of our system. With thermal pumping schemes 
however, a genuine steady state is achieved only for an average over many trajectories. 
For a given set of parameters, each trajectory differs not only in the actual relative phase 
established between the two traps, but more importantly, in the instantaneous atom numbers 
as a function of time. As the trajectory simulation proceeds, the occupancy of each trap 
exhibits thermal fluctuations which lead to time variations in other properties of the system 
such as the conditional visibility. Only the time-averaged properties approach a true steady 
state. As discussed in section 0, we treat two cases: two-way pumping in which atoms may 
be exchanged between the bath and trap in both directions, and one-way pumping in which 
atoms can only move from the bath into the trap. We begin with representative plots of 
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the visibility as a function of time for a single trajectory with no collisions (k = 0) . The 
visbility for two-way pumping is shown in Fig. There are initially n = 100 atoms in each 
trap and the pumping rate is chosen to balance the detection rate. On average, 2n atoms 
are detected in a time 7^ = 1. The visibility is extremely noisy with frequent fluctuations of 
order 1. Our simulations show the occupancies of the traps also display large fluctuations 
as would be expected for coupling to thermal baths In particular, a zero in the visibility 
is always associated with a zero in one or other of the atom numbers. The visibility for 
a typical trajectory with one-way pumping is shown in Fig. ^d. In this case, there are 
again large fluctuations but on a much longer time scale. This difference has a simple origin 
mentioned earlier: for the one-way case, on average one atom is added to the system for 
each atom detected and so all the atoms are replaced once in time 'jt = 1. In the two-way 
case, n atoms are exchanged with the baths for each atom detected, and so in 7^ = 1, all 
the atoms are replaced n times over and we expect a correspondingly shorter time scale for 
the fluctuations. As a comparison, in Fig. 0c we show a trajectory for a "regular" pumping 
model in which atoms are dripped into the trap at a constant rate to replace those lost by 
detection. In this case, the collisional rate is k = O.57, but the visibility shows a much 
improved response than in the (coUisionless) thermal pumping cases, indicating the severe 
influence of the thermal pumping. 

In fact, the visibility is degraded by the number fluctuations in two distinct ways. When 
the occupancy of one of the traps falls due to a fluctuation to within a few times (T„ of zero, 
the extreme terms in the entangled state are removed, the number distribution narrows and 
the visibility falls. In particular, if one of the traps is completely emptied (as occurs several 
times in Figs. 0a and b,) the state is then a pure number difference state and any relative 
phase is completely destroyed. The visibility can of course be restored once fluctuations 
increase the atom number again, but there is no relation between the new relative phase 
and the phase before the trap was emptied. 

Even when both traps have (n) ^ the visibility is reduced according to the number 
difference between the traps. This is obvious — if one trap has significantly more atoms than 
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the other, then we can predict with better than 50 % accuracy from which trap the next 
atom will be detected, and the visibility must fall accordingly. We can calculate this effect 
simply as follows. In our system, the atom numbers experience thermal fluctuations in time 
due to the pumping. Suppose for a moment we have a different situation in which there is 
no pumping, and we perform a series of detection runs with a thermal distribution in the 
initial trap numbers and measure the visibility after a well-defined relative phase has been 
set up (but before the traps are significantly depleted). If we picture the condensates as 
coherent states with some relative phase: 

I )i = \/^exp (i0i) , I )2 = v^exp (i02) , (15) 

the expected visibility is just /3 = ^^fn^n^ {n\ + ^2) which is a familiar expression for 
optical fields. Defining the relative occupancy 

/ = (16) 

we have 



/?(/) = Vl-/^. (17) 

More correctly, we should derive Eq. (|l^ directly from the entangled state description of 
the twin-trap system. In general, the fringe visibility is given by 



/5=k^'^|VW^, (18) 

where U*-^-*! is the normalized correlation function pT|. While in general, 1 5'*-^'' I is 



not easily evaluated for the entangled states with which we are concerned, it can 



be shown for example that for the projected two-mode coherent state = 
SfcLo Oi^(^^~^ / \/ k\ {N — A;)! I A;, N — k) , which is the most natural expression of a state with 
relative phase with fixed total atom number N, \g^^^\ tends to unity in the limit of large N, 
Figures |^ and b test Eq. (|T^) in the form of scatter plots of the points (/ (t) ,P (t)) 
for two-way and one-way pumping respectively with the same parameters as Figs ^ and 
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b. The prediction ([T7| ) is indicated by the black squares in each plot. The correlation is 
clearly much stronger in the one-way case. This difference is entirely due to the difference 
in time scales discussed above. If we are to have a well-defined phase, the number must be 
partially uncertain. Indeed, we see later that in the presence of collisions the average state 
has moderate number-squeezing but with a variance of the same order as a coherent state. 
So for a reasonable visibility we should require an entangled state of order a = O (2y^) 
terms. Such an entangled state is set up by the same number of detections and requires 
a time of order = l/v^7- Now for two-way pumping, the time scale for replacement of 
all the atoms once over is Tj. = 1/727 ^ ^e- Hence, the exchange of atoms with the baths 
occurs faster than an entangled state of a particular phase can be constructed and we may 
expect a reduced visibility. There can be only a weak correlation between the instantaneous 
visibility and the instantaneous relative occupancy / , and the visibility is generally lower 
than the optimum given by Eq. (p!?!). With one-way pumping however, the time scale for 
replacement of all the atoms is larger by a factor n. The visibility is able to keep up with 
the drift in number and is then limited only by Eq. ([T7|). 

We can also calculate the mean visibility over time /3, for an arbitrary pair of mean atom 
numbers {rii) and (722) . Again, we think of an ensemble of runs with no pumping and a 
thermal distribution of initial states and \n2) ■ The mean visibility over many runs is 
the weighted average of (3 (/) over the probability distribution [see Eq. ([T6|) ] 



where Pn. (rii) = — log {•ji) 7^ ' are the probability distributions (in the continuous limit) of 
atom number for thermal distributions with mean number and 7^ = n^/ (n^ -|- 1) . For the 
case where the mean numbers are the same, P (/) is uniform and (3 = tt/A [^]. Otherwise 
we find 




(19) 



27rlog (71) log (72) 
[log (71/72)]^ 




1 



1 



(20) 
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which for ni, n2 S> 1, gives 

(1 + v^) 

with p = ni/n2- In the pumped twin-trap setting, we also have thermal distributions in the 
atom number which occur not from run to run, but over time in a single trajectory, so it is 
reasonable to hope that the above analysis may still apply. 

In Fig. ^ we show the average visibility for a thermal distribution as a function of the 
mean atom number ratio p given by Eq. (^if) . The plotted points show the time-averaged 
visibility calculated from trajectory simulations with one-way pumping and ni + n2 = 200. 
Error bars are shown at 1 standard deviation. As expected, the mean visibility falls with 
increasing disparity in the mean atom number. 



A. One-way pumping and output couplers 

We have seen that the one-way pumping process shows significantly higher visibility 
than the two-way pumping. From this point on, we restrict our attention to the one-way 
model and add the effect of an output coupler from each trap. In our results we find two 
distinct regimes according to the length of the simulations. Figure |^ shows the visibility as 
a function of time averaged over 200 trajectories for a simulation with k = 0, z/j = and 
an initial state |ni, ^2) = |100, 100) . The one-way pumping rate was chosen to maintain the 
mean population at n = 100 in each trap. There are clearly two regimes: for '-ft ^ n, the 
mean visibility shows a steady decline, while for 'jt n, the visibility tends to a steady-state 
value of 7r/4 consistent with the calculation of the previous section. In the initial stage of a 
particular run, the populations of the traps become decorrelated due to the thermal nature 
of the pumping until they are completely uncorrelated and the time-averaged visibility is 
7r/4. The time for this decorrelation varies from run to run, having a characteristic length of 
'jt ^ n. Thus the average over many trajectories shows a gradual decline until all members 
of the ensemble are likely to be decorrelated. We are thus led to examine the behavior of the 
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system in the two regimes •yt <^ n, when the trap populations are hkely to be quite close, 
and '-ft ^ n, when there is no correlation between the populations. We treat these two cases 
in turn. In all cases, we start our simulations with the initial state \ni,n2) = |100, 100) and 
calculate quantities averaged over 200 trajectories. 

• 'jt <^ n : For the short-time regime, we arbitrarily choose 7t = 4 to show results. Fig- 
ure l^a shows the visibility as a function of k again averaged over five trajectories. As 
expected, the visibility decreases with increasing collisions which increasingly disrupt 
the relative phase 0]. We have also performed simulations with a range of output 
couplings from i/j = to z/j = 7. This leads to a small decrease in (3 (of less than 0.025 
for the strongest coupling). This effect is simply a result of the fact that the pumping 
rate is increased to balance the additional loss of atoms and so the trap populations 
decorrelate faster. The nature of the average state of the system is indicated in Fig. ^b. 
Shown are the widths of the number distribution cr„ (dotted line) and phase distribu- 
tion (dot-dashed) and p: the root mean square of the product of the two (solid). 
The filled circles denote the actual simulations performed. The number distribution 
clearly narrows strongly with increasing coUisional rate while the phase distribution 
spreads as collisions degrade the relative phase. The simulations with non-zero output 
coupling (not shown in Fig. Po) produced a reduction of less than 5 % in the number 
width and no discernible change in the other parameters. Note that for zero collisions, 
the product of the widths (solid) is unity indicating a minimum uncertainty state. 
Further, for all values of k, the number width cr„ < = 10, which is the width we 
would expect if the state was a projection of a coherent state onto a basis of fixed 
total atom number. The real state is thus quite strongly number squeezed. This in 



consistent with recent analytic work by Dunningham et al 1^5^11 using a Bose-broken 
symmetry model. They show that in the limit of a large collisional rate, the true state 
of the condensate is the amplitude-squeezed state that minimizes number fluctuations. 

7t S> n : In the large time regime, the mean visibility has no dependence on the output 
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coupling rate — once the atom numbers are completely uncorrelated, the precise rate 
at which atoms enter the trap is irrelevant. The number and phase widths shown 
in Figure |^ show very similar trends to the short time case. Note that even with the 
uncorrelated trap numbers, the state is still minimum uncertainty for k, = 0, indicating 
that the pumping rate is low enough for the visibility to adjust to changes in number. 
Again, other simulations showed that the only effect of output coupling was to reduce 
the number width by a few percent. 



V. APPLICATION TO COLLAPSES AND REVIVALS 

In this final section, we consider the application of pumping processes to the interesting 
phenomenon of collapses and revivals in the relative phase. Several authors have shown that 
if a relative phase is prepared by detection and the entangled state subsequently evolves 
purely under the influence of the interatomic collisions, the visibility of the phase experiences 
recurrent collapses and revivals of period tt/k due to the differential rate of phase rotation 



in the entangled state [^,§,0,|ri|] . A demonstration of collapses and revivals of the phase, 
perhaps through light scattering experiments, would be a significant result in BEC physics. 
It is interesting to consider how the collapses and revivals are affected by pumping and 
leaking of atoms through the traps. Naively, we might expect that the oscillations would be 
destroyed by the time all the atoms had been replaced a few times over. In fact, we have 
found the collapses and revivals to be remarkably robust to pumping processes. 

In Fig. I^a, we show the visibility for a single trajectory without pumping or output 
coupling in which there are initially 200 detections from a total of 1000 atoms, followed by 
a period in which the system evolves only under the influence of collisions with n = 0.25. 
The oscillations in the visibility are clear. Fig. ^3 shows a trajectory with the same number 
of detections and collision strength, but with a continual flushing of the trapped atoms by 
pumping and output coupling. On average, all the atoms are replaced in a period 7t = 1 and 
the atom numbers exhibit large fluctuations (Fig. |c). Despite this, the collapse and revivals 
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persist for a considerable period and only disappear when the atom number in trap 2 (thick 
line in Fig. ^) approaches zero at 'jt ^ 250. If the trajectory is such that neither atom 
number approaches zero, the oscillations may continue much longer still. As the simulation 
progresses however, while the period of the revivals is unchanged, the peaks broaden — the 
collapse time increases. This is associated with the gradual reduction in the width of the 
number difference of the initial entangled state due to the repeated addition and removal of 
atoms indicated in Fig. |^d. Essentially, every addition or removal of an atom through an 
output coupler tends to drive the state to a narrower number distribution. In their treatment 
of collapses and revivals for a fixed number of atoms, Wong et al. § have shown that the 
visibility during a collapse should decay according to 

V oc exp {-2(j\kH'^) , (22) 

where a a is the width of a Gaussian approximation to the coefficients 

A{k) = |cfcCfc_i| ^/{n — k + 1) {n — m + k), (23) 

and there have been m detections from an initial state of n atoms in each trap. For a broad 
distribution, and m -C n, to lowest order we have A{k) oc |c^| , so that in our notation 
aA ~ o'n/2. The black squares in Fig. ^d are estimates of cx^ calculated from the collapse 
widths in Fig. §3 using Eq. (|22D . The agreement with the directly measured values for the 
width of the number distribution (solid line) confirms that the increase in the collapse time 
is due purely to the change in (T„. 

Figure §3 also shows a variation in the height of the visibility peaks. Note that the 
variation is not monotonic, an effect we have found to be generally true. A natural guess is 
that the peak heights are associated with the relative occupancy / = (rii — ^2) / {rii + 712) , 
which we found led to a maximum visibility for systems where the detections are not stopped 
in section |V[ We have tested this using scatter graphs of the peak visibility similar to those 
in Fig. 1^. We find a moderate confirmation of the connection. In cases for which the minima 
of the visibility remain small, there is a strong correlation between the peak visibility and 
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the quantity /. In other cases, such as that in Fig. §3 for 7^ > 150, for which the minima 
are significantly greater than zero, the correlation is poor and we conclude that the pumping 
process has produced an additional degradation of the state beyond that implied simply by 
the mean number difference. 

VI. CONCLUSION 

In this paper, we have studied the steady-state behavior for two pumping scenarios to 
show how an ongoing measurement process can generate a phase coherence between atoms 
derived from thermal baths, even in the presence of phase diffusion due to atomic collisions 
within the traps. We find important qualitative differences between systems with two-way 
and one-way pumping, the phase coherence being substantially improved for the one-way 
case. Systems displaying collapses and revivals of the condensate phase should provide an 
opportunity for examining the time- dependent effects induced by pumping. We remark 
finally that a natural extension to our model would be the inclusion of extra trap levels that 
would allow for line narrowing and genuine laser action. 
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FIG. 1. Geometry of pumped twin-trap system. The straight sohd arrows indicate the detection 
at rate 7; the dotted arrows, the exchange of atoms with the reservoirs; the curved arrows, output 
couphng of the trapped atoms. 

FIG. 2. a)-b) Visibility /? as a function of time for a thermally pumped system with mean 
occupancy n = 100 for each trap and k = 0. a) Two-way pumping, b) one-way pumping, c) 
Visibility for a regularly pumped model with k = 1.0. 

FIG. 3. Scatter plot of /? as a function of relative occupancy / for a) two-way, and b) one-way 
pumping. There are 10000 points shown. Black squares indicate the relation Eq. ( p!7| ) 

FIG. 4. Mean visibility /3 as a function of the atom number ratio p. Error bars indicate 1 
standard deviation errors in time-averaged simulations. 

FIG. 5. Visibility for one-way thermal pumping with no collisions and no output coupling. The 
mean atom number in each trap is 100. 

FIG. 6. Averaged state parameters as a function of collision rate for one-way pumping in 
short-time regime, a) Visibility and b) (T„ (dashed), a^p (dot-dash), and p (solid). 

FIG. 7. Averaged state parameters as a function of collision rate for one-way pumping in 
long-time regime: cj„ (dashed), a^p (dot-dash), and p (solid). 
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FIG. 8. a) Visibility for collapses and revivals of relative phase with no pumping or output. 
Initially 200 detections were made from 1000 atoms, b) Visibility, c) atom numbers and d) an for 
the same parameters with pumping and output rates such that all atoms are replaced on average 
once in a period 'yt = 1. 
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